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We study the process of spinodal decomposition in a scalar quantum field theory 
that is quenched from an equilibrium disordered initial state at T% > Tf to a final 
state at Tf ~ 0. The process of formation and growth of correlated domains is 
studied in a Hartree approximation. We find an approximate scaling law for the size 
of the domains rj v^£o at long times for weakly coupled theories, with £o the 

zero temperature correlation length. 
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The process of phase separation through spinodal decomposition is well understood 
within the context of classical non-equilibrium statistical mechanics [1|-|3]]. When the quench 
is at criticality, it is primarily associated with the onset and growth of unstable long- wave- 
length modes (we will not discuss nucleation here). 

Spinodal decomposition is conjectured to have taken place during phase transitions in the 
early universe typically described by scalar quantum field theories Thus it becomes 

an important issue to understand and describe the mechanism of phase separation in quan- 
tum field theory with the definite motivation provided by cosmological phase transitions. 
An attempt to study some of these issues was reported by Mazenko and collaborators J7j. 
Understanding the process of phase separation in quantum many-body systems may also 
prove relevant within the context of mesoscopic systems and quantum phase transitions in 
condensed matter. In this article, we introduce the techniques of non-equilibrium quantum 
statistical mechanics to study the process of phase separation in a typical scalar field theory. 

From the outset we recognize fundamental differences between classical theories of spin- 
odal decomposition, and the description of phase separation in a quantum system. In the 
former, the equations of motion are purely dissipative, whereas in the latter, the Heisenberg 
operator equations of motion are second order and thus time reversal invariant. Further- 
more, thermal fluctuations are incorporated in the classical description via a Langevin noise 
term, typically uncorrelated and satisfying the fluctuation-dissipation relation, whereas in 
the quantum description, both thermal and quantum fluctuations are present in the initial 
density matrix that describes the initial ensemble. 

We model the physical situation of a "critical quench" via a time dependent Hamiltonian 

Hit) = Jj 3 x |in 2 (x) + i(V$(a;)) 2 + \m>(t)<S>\x) + ^ 4 (*)} (0-1) 
m 2 (t) = m 2 6(-t) - m}e{t) (0.2) 

™? = M 2 [3!-1] ! rn) = n 2 [±-%] (0-3) 

c c 

with fi 2 > ; Ti > T c ; Tf ~ 0. The initial state of the system (at t < 0) is assumed to be 
described by an equilibrium density matrix at the initial temperature Tj. 
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Pi = 

Hi = H(t < 0) 

In the Schroedinger picture, the density matrix evolves in time as 

p(t) = U(t)p t U~\t) 



(0.4) 
(0.5) 



(0.6) 



with U(t) the time evolution operator. In the present case, the order parameter 
(Trp(t) f n d 3 x§(x)) obeys a Heisenberg equation of motion and is not conserved. The ex- 
pectation value of operators may be computed by introducing the generating functional 



Z[J+, J-} = TrUiT - i/3i, T)U(T, T'\ J~)U(T', T; J + " 



(0.7) 



with T — > — oo; ; T' — > oo. It is found 



T' 



Z[J + , J-] = exp U / dt C int (-i5/5J+) - £ mt (i5/5J-) 



x 



with G a b the Green's function on a contour M . The quantities of interest are 



(0.8) 



S(r;t) = mr,t)<S>(0,t)) 
d 3 k 



S(f; t) 



(2tt)' 



(0.9) 
(0.10) 



The equal time correlation function at zeroth order ( "tree level" ) is thus found to be 
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<*s(t)*_s(t)> 

with the mode functions obeying 

d?_ 
dt 2 



2cj<(A;) 



W+(t)W fc -(t)coth[Au;<(fc)/2] 



(0.11) 



+ k 2 + m 2 (t) 



= 



(0.12) 



with m 2 (t) given by ( |0.2| ). 

The boundary conditions on the homogeneous solutions are 



U£{t < 0) = e 



=pioj< (k)t 



k 2 + m 2 



(0.13) 
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corresponding to positive frequency (particles) and negative frequency (antiparticles) 
(U£(t);Uk(t), respectively). 

The "free field" mode functions are easily found. For t > 0, they consist of stable modes 
(k 2 > mV) and unstable modes (k 2 < to 2 ). These unstable modes are responsible for the 
growth of correlations. The zeroth order equal time correlation function becomes 

= coth[Ao,<(A;)/2] (0.14) 

for t < 0, and 

= {[1 + ^ k B k [cosh{2W{k)t) - 1]] 0(m 2 - k 2 ) 

+ [1 + 2a k b k [cos(2cu > (k)t) - 1]] Q(k 2 - m))} coth[/3^<(A;)/2] (0.15) 



with u > (k) = \Jk 2 — mj and W(k) = \Jm 2 — k 2 , for t > 0. 

The first term, the contribution of the unstable modes, reflects the growth of correlations 
because of the instabilities and will be the dominant term at long times in this approximation. 

It is convenient to introduce the dimensionless quantities 

k r2 m 2 \Tf-T 2 } _ , , 

k= — ; L 2 = — j = f^ c -\ ■ r = m f t ; x = m f f (0.16) 



rn/ mj T 2 — TJ 

and the critical temperature T 2 = 24/i 2 / A in terms of which the "tree level" subtracted 
structure factor S^°'(k,t) — S^(k, 0) = (1/to/)«S^(k, t) becomes 

S i0) ^r) = | | (p) (l + ^) [cosh(2^r) - 1] (0.17) 




lu 2 = k, + L ; W K = 1 — K (0.18) 



To obtain a better idea of the growth of correlations, it is convenient to introduce the scaled 
correlation function 

\ A f m f k 2 dk sm(kr) t „,, , 

The reason for this is that the minimum of the tree level potential occurs at A$ 2 /6m 2 = 1, 
and the inflexion (spinodal) point, at A$ 2 /2m 2 = 1, so that X>(0, r) measures the excursion 
of the fluctuations to the classical spinodal and beyond as the correlations grow in time. 



At large r (time), k, 2 S(k,, t) has a sharp peak at k s = 1/y/r with amplitude exp[2r]/r 



Restoring dimensions, and recalling that the zero temperature correlation length is £ (0) = 
l/v2//, we find that for Tf ~ the amplitude of the fluctuation inside a "domain" ($ 2 (t)), 
and the "size" of a domain £d(£) grow as 



The presence of the instabilities precludes a well-defined perturbative expansion. Con- 
sider a one loop contribution to the equal time correlation function. The "external legs" 
obtain a contribution from the unstable modes, but in the loop integral, the integration 
over the momenta also includes a contribution from the unstable modes. It is clear that 
eventually the one-loop correction dominates and perturbation theory breaks down, even 
for the case of very weak couplings. This feature will persist to all orders in a perturbative 
expansion. The dynamics of the phase transition cannot be studied in perturbation theory. 

Our non-perturbative approach is based on a Hartree approximation, which is similar to 
the early approach of Langer |J for classical spinodal decomposition. It is implemented by 
the replacement 



(see figure 1). We find for x < r and Tf « 




(0.20) 



(0.21) 




(0.22) 



m 2 (t) 



m 2 (t) + -($ 2 (t)) 



(where we used spatial translational invariance). 



This leads to the self consistent set of equations 



dt 2 



+ k 2 + m 2 (t) + ^ 2 (t)) Ut 







(0.23) 




(0.24) 
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The composite operator ($ 2 (r, t)) needs one subtraction and multiplicative renormal- 
ization. The subtraction is absorbed in a renormalization of the bare mass, and the mul- 
tiplicative renormalization into a renormalization of the coupling constant. The Hartree 
approximation provides a self-consistent non- perturbative scheme that sums an infinite 
series of Feynman diagrams [O . For t < there is a self-consistent solution given by 



($ 2 (t)> = ($ 2 (0)> ; Ut(t)=exp[^u; < (k)t} 
^(fc) = P + m 2 + ^ + ($ 2 (0)) = P + m 2 



(0.25) 



(0.26) 



and mf R = /^[(Zf/T 2 ) — 1]. For t > we subtract the composite operator at t = absorbing 
the subtraction into a renormalization of m 2 which we now parametrize as m 2 R = fi R [l — 
(T 2 /T 2 )]. This choice of parametrization only represents a choice of the bare parameters. 
The logarithmic multiplicative divergence of the composite operator will be absorbed in a 
coupling constant renormalization consistent with the Hartree approximation [|l(J . However, 
for the purpose of understanding the dynamics of growth of instabilities associated with 
the long-wavelength fluctuations, we will not need to specify this procedure. After this 
renormalization, the Hartree equations read 

d 3 k 1 



imt)} - mo)}] 



(27r) 3 2cu<(A;) 



[U£(t)Ut(t) - 1] coth[A^<(fc)/2] 



^ + fc 2 + m 2 ,(t) + ^(($ 2 (t))-($ 2 (0 



u£(t) = o 



(0.27) 



(0.28) 



J^2 



- 1 



©H)-/4 



T 2 

zL 

J^2 



e(t) 



(0.29) 



with Tj > T c and Tf <ti T c . With the self-consistent solution and boundary condition for 
t < 



[($ 2 (t < 0)) - ($ 2 (0))] = ; Ut{t < 0) = exp[^<{k)t] 
u K (k) = \Jk 2 + m% 



iR 



(0.30) 
(0.31) 
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This set of Hartree equations is extremely complicated to be solved exactly. However it 
accounts for the process of coarsening Consider the equations for t > 0, at very early 
times, when (the renormalized) ($ 2 (t)) — ($ 2 (0)) ~ the mode functions are the same as in 
the zeroth order approximation, and the unstable modes grow exponentially. By computing 



the expression (|0.27|) self-consistently with these zero-order unstable modes, we see that the 



fluctuation operator begins to grow exponentially. 

As ((& 2 (t)) — ($ 2 (0))) grows larger, its contribution to the Hartree equation tends to 
balance the negative mass term, thus weakening the unstabilities, so that only longer wave- 
lengths can become unstable. Even for very weak coupling constants, the exponentially 
growing modes make the Hartree term in the equation of motion for the mode functions 
become large and compensate for the negative mass term. Thus when 



2m },R 



$ 2 (t)) - ($< 



the instabilities shut-off, this equality determines the "spinodal time" t s . The modes will 
still continue to grow further after this point because the time derivatives are fairly (expo- 
nentially) large, but eventually the growth will slow-down when fluctuations sample deep 
inside the stable region. 

After the subtraction, and multiplicative renormalization (absorbed in a coupling con- 
stant renormalization), the composite operator is finite. The stable mode functions will 
make a perturbative contribution to the fluctuation which will be always bounded in time. 
The most important contribution will be that of the unstable modes. These will grow expo- 
nentially at early times and their effect will dominate the dynamics of growth and formation 
of correlated domains. The full set of Hartree equations is extremely difficult to solve, even 
numerically, so we will restrict ourselves to account only for the unstable modes. From the 
above discussion it should be clear that these are the only relevant modes for the dynam- 
ics of formation and growth of domains, whereas the stable modes, will always contribute 
perturbatively for weak coupling after renormalization. 

Introducing the dimensionless ratios (|0.16| ) in terms of mf t R ; m^, (all momenta are now 



expressed in units of mf :R ), dividing ( |0.28|) by rn 2 Rl using the high temperature approx- 
imation coth[/3iUj < (k) /2] « 2T i /uj < (k) for the unstable modes, and expressing the critical 
temperature as T 2 = 24/i|,/ X R , the set of Hartree equations (|0.27| , |0.28| ) become the follow- 
ing integro-differential equation for the mode functions for t > 



dr 2 



+ q 2 - 1 + / dp 



P 



p 2 + L 2 



K{r)U-{r) 



(0.32) 



with the boundary conditions (p.30|) for t < and 

T; 
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(0.33) 



4vr 2 [T 2 -T 2 ]5 

The effective coupling (p.33|) reflects the enhancement of quantum fluctuations by 
high temperature effects; for Tf/T c « 0, and for couplings as weak as X R ~ 10~ 12 , 
g 10 _7 (Tj/T c ). This value of the coupling has particular significance in inflationary mod- 
els and arises from bounds on density fluctuations £§-|6|. The equations (p.32|) may now be 
integrated numerically for the mode functions; once we find these, we can then compute 
the contribution of the unstable modes to the subtracted correlation function equivalent to 



V^ hf \x,t) 
3V^ hf \x,t) 



6m 2 R 



($(f,t)$(0,t))-($(r,0)$(0,0)) 



g dp 



P 



\ sm(px) 
j) 2 + L 2 J (px) 



u;{t)u;{t)-i 



(0.34) 
(0.35) 



In figure (2) we show 3(T> HF (0,r) — V HF (0,0)) (solid line) and also for comparison, its 
zeroth-order counterpart 3(£>(°>(0,t) - P^(0,0)) (dashed line) for X R = KT 12 , T;/T c = 2. 
(This value of the initial temperature does not have any particular physical significance 
and was chosen as a representative). We clearly see what we expected; whereas the zeroth 
order correlation grows indefinitely, the Hartree correlation function is bounded in time and 
oscillatory. At r w 10.52 , 3(V^ HF \0,r) - V^ HF \0,r)) = 1, fluctuations are sampling field 
configurations near the classical spinodal, fluctuations continue to grow, however, because 
the derivatives are still fairly large. After this time, the modes begin to probe the stable 



region in which there is no exponential growth. At this point 



2ni 



f,R 



-(($ 2 (r)> - $ 2 (0)», 



becomes small again because of the small coupling g pa 10~ 7 , and the correction term 
becomes small. When it becomes smaller than one, the instabilities set in again, the unstable 
modes begin to grow and the process repeats. This gives rise to an oscillatory behavior 
around 2r ^f ((<l )2 (r)) — $ 2 (0)) = 1 as shown in figure (2). We clearly see that for very weakly 
coupled theories, the zeroth order correlation function provides a fairly good approximation 
to the Hartree correlations up to the "spinodal time" . Thus for very weakly coupled theories 
correlation functions will be approximately given by ( p.20| , |0.21| ) and this permits us to find 
an approximate result for the spinodal time (at which fluctuations begin probing the stable 
region) . 



t. . 
r, = — = « — In 

v^(0) 



'3A' 

471-3 



/ 



Tj \3 
^2T C > 

[R-i] 



(0.36) 



It is remarkable that the domain size scales as £d(£) ~ £2 just like in classical theories 
of spinodal decomposition, in which the order parameter is not conserved, as is the case 
in this relativistic scalar field theory, but certainly for completely different reasons. At the 
tree level, we can identify this scaling behavior as arising from the relativistic dispersion 
relation, a situation very different from the classical description of the Allen- Cahn- Lifshitz 
fll| theory of spinodal decomposition based on a Time- dependent Landau Ginzburg model. 
For strong coupling, the Hartree result and the zeroth-order result depart at very early times 
fl"2|j . It is well known within the context of classical spinodal decomposition that the Hartree 
approximation is not correct at intermediate and long times. We are currently studying a 
consistent treatment beyond the Hartree approximation. Details of the calculation and 
possible extensions will be presented elsewhere |JL2|| . Of particular importance will be the 
study of the interface dynamics, known to be the relevant description in classical theories 
0. 
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Figure Captions: 
Figure 1: 

(rexp[-2r])K5(K,r) for A = 10~ 12 ; T;/T c = 2 and r = 10, 13, 16, 19. 
Figure 2: 

Hartree (solid line) and zero order (dashed line) results for ^r({^ 2 (j)) — ($ 2 (0))) = 
3P(0,r), for A = 10~ 12 , g = 2. 
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